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ABSTRACT 

The well-known relationship between metallicity and luminosity in active galactic nuclei (AGNs) is addressed 
by introducing new metallicity measurements (based on the method of Hamann & Ferland, hereafter HE) for a 
sample of narrow-line Seyfert 1 galaxies (NLSls). Our new results, based on a sample of 162 AGNs, including 
nine NLSls, indicate that while broad-line AGNs trace a metallicity-luminosity power law with an index of 
~ 0.2, NLSls deviate significantly from this relationship at low luminosities. Adopting the HE method based on 
the N V/C IV line ratio, we find that NLS 1 metallicities are similar to those of some high-redshift, high-luminosity 
quasars. We also examined the N IV]/C IV line ratio and compared it with N V/C IV in a sample of 30 sources 
including several NLSls. We find that the two do not give a consistent answer regarding the N/C abundance 
ratio. This result is marginal because of the quality of the data. We suggest two alternative explanations to these 
results: 1) The HE metallicity-luminosity dependence is not a simple two-parameter dependence and there is an 
additional hidden variable in this relationship that has not yet been discovered. The additional parameter may 
be the accretion rate, the age of the central stellar cluster or, perhaps, something else. 2) The strong line ratios 
involving N V A 1240 suggested by HE are not adequate metallicity indicators for NLSls and perhaps also other 
AGNs for reasons that are not yet fully understood. 

Subject headings: galaxies: abundances — galaxies: active — galaxies: nuclei — galaxies: Seyfert 



1. THE METALLICITY-LUMINOSITY RELATIONSHIP IN 
ACTIVE GALACTIC NUCLEI 

Studies of emission lines in active galactic nuclei (AGNs) 
indicate that metallicities are typically near the solar value in 
their broad- line region (BLR). Accurate determinations are dif- 
ficult to obtain since the line ratios depend on unknown densi- 
ties and optical depths (Netzer 1990). However, Shields (1976) 
showed that some BLR abundances can be derived from the ob- 
served line ratios of weak nitrogen, carbon, and oxygen lines, 
almost independent of the physical properties of the gas (see 
Hamann et al. 2002 for more details). One such line ratio is 
N IV] A1486/C IV A1549 with a theoretical value of - 0.045 
for solar metallicity. Others are N III] A1750/O III] A1663, 
O III] A1663/C IV A 1549 and N III] A1750/C III] A 1909, which 
can determine the N/0, O/C, and N/C abundances, respectively. 
This method has been used in several AGN studies (e.g., Bald- 
win & Netzer 1978; Osmer 1980; Uomoto 1984), but there are 
practical limitations due to the weakness of these lines. 

An important development in this area was achieved by 
Hamann & Ferland (1993, hereafter HF93; see also Hamann 
& Ferland 1999), who suggested alternative abundance indi- 
cators that are somewhat model dependent but much easier to 
obtain observationally. In particular, they have shown that the 
N V A1240/C IV A 1549 and N V A1240/He II A 1640 line ra- 
tios represent the overall BLR metallicity. They also claimed 
that BLR metallicity tends to grow with AGN luminosity up 
to ^ IOZ0, thus implying a metallicity-luminosity (hereafter 
Z-L) relationship, in analogy with the mass-metallicity rela- 
tionship observed for elliptical galaxies. At one extreme end 
of this relationship one finds the luminous and high- redshift 
quasars as AGNs having the highest metallicities. In several 
cases, those quasars display relatively narrow UV emission 
lines (e.g., Osmer 1980; Warner et al. 2002). 

The low-luminosity regime in the HF93 diagram is also oc- 
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cupied by narrow-line Seyfert 1 galaxies (NLSls) that have 
been mostly left out from previous abundance analyses. These 
NLSls are defined by their extremely narrow optical per- 
mitted emission lines (FWHM ^ 2000 km s~') in comparison 
with normal broad-line AGNs (BLAGNs; Osterbrock & Pogge 
1985). NLSls show extreme AGN properties; their optical 
emission lines put them at one extreme end of the Boroson 
& Green (1992) primary eigenvector, and they tend to display 
unusual behavior in other wave bands, especially in the X-ray 
(e.g.. Boiler, Brandt, & Fink 1996; Leighly 1999a, 1999b). A 
possible explanation for the peculiar properties of NLS 1 s is that 
they have relatively low black hole (BH) masses for their lumi- 
nosities and hence a very large L/LEdd- Since only a handful 
of NLSls had their BH mass measured directly using rever- 
beration mapping techniques (Peterson et al. 2000), this is not 
yet fully confirmed. Another extreme NLS 1 property was re- 
cently suggested, namely, unusually high metallicities (Mathur 
2000 and references therein). However, to date, this evidence 
remains scarce, and no systematic abundance study has yet been 
carried out. 

In this study we present for the first time metallicity (a la 
HF93) measurements for a sample of NLSls. In § 2 we de- 
fine the sample properties and data analysis. In § 3 we describe 
our new results on the Z-L relationship in AGNs and attempt 
to answer two related questions: (1) Do NLSls have higher 
metallicities compared with BLAGNs for a given luminosity? 
(2) Are the higher metallicities, assumed for NLSls related to 
the fundamental physical properties that drive the NLS 1 phe- 
nomenon? 



2. SAMPLE PROPERTIES AND DATA ANALYSIS 

We selected 162 type 1 AGNs, including nine extreme 
NLSls, that had either published optical or UV (HST) HF93 
line ratios or an UV (HST-archived) good-quality spectrum 
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Fig. 1.— Two typical line-fitting examples: A BLAGN PG 1202+281 (top 
panels) and an NLSl KUG 1031+398 {bottom panels). In each continuum- 
subtracted spectrum, the Lya region {left panels) and the C IV region {right 
panels) were fitted as explained in the text. In each panel a thin solid line rep- 
resents the spectrum; dotted lines represent the individual Gaussian profiles 
that are then summed to form the entire fit to the spectrum {thick solid lines). 

covering the rest-frame range of ^ AA1200-1700. We also in- 
spected \JV(1\JE) spectra of 1 1 extreme NLSls, but tiie poor 
signal-to-noise ratios resulted in line flux uncertainties that 
were too large to be included in our study (although the general 
trend was similar to the nine NLSl spectra taken from HST). 
A redshift of z « 1 .6 defines the transition from low-z objects, 
for which the relevant emission lines are detected in the UV, 
to those objects for which N V is detected in the optical band 
and are hereafter referred to as high-z quasars. We consider 
NLS 1 s as AGNs having FWHM(H/3) < 1 500 km s-Un order to 
concentrate on the more extreme objects of this class ^. 

Among the 162 objects, 67 have published optical spectral 
data and are regarded as high-z and unknown FWHM(H/3) ob- 
jects (Osmer & Smith 1976, 1977; Baldwin & Netzer 1978; Os- 
mer 1980; Uomoto 1984; Baldwin, Wampler, & Gaskell 1989; 
Baldwin et al. 1996; Dietrich & Wilhehn-Erkens 2000). All 
those include detailed line intensity measurements. The rest, 
95 low-z objects, including the nine extreme NLSls, have UV 
(HST) spectra. Of these, 44 have published line intensities (Re- 
ichert et al. 1994; Laor et al. 1994, 1995, 1997; Wills et al. 
1995). Emission-line fluxes of 111 objects were thus readily 
available for our analysis. For the remaining 51 objects, we 
obtained the HF93 line ratios by performing line and contin- 
uum measurements on the archived spectra by applying a multi- 
Gaussian component fit using the NGAUSSFIT task in IRAF"*. 
For Lya and C IV A1549, two Gaussian components were fit- 
ted, while for N V A1240, He n A1640, and O m] A1663 we 
fitted a single Gaussian. Two typical fits of this kind appear 
in Figure 1 . The individual fits produced uncertainties on the 
line fluxes that were much smaller compared with the uncer- 
tainties introduced by other unknowns, such as the precise con- 
tinuum placement and slope, line contamination, etc. In partic- 
ular, the Si II A 1264 complex might be expected to contribute 
at most ~ 20% (Laor et al. 1997) to the N V A1240 intensity 
only in AGNs with FWHM(N V) > 5000 km s-\ In Ught of 
the underestimated uncertainties, we compared published flux 
values with new measurements that we performed for several 




10" 10'" 
vL.@1450A(L3„J 

Fig. 2. — HF93 metaUicity indicators, Nv/Civ (top) and Nv/Hell (bot- 
tom), as a function of luminosity. Empty circles mark BLAGNs, and filled 
circles marks NLSls; solid lines represent the BLAGN best-fit Z-L slope; 
large squares with error bars represent average line ratios in bins of 0.5 in log 
vLi, of the entire data set. Note the significant deviation of the low-luminosity 
bins from the straight line (top) owing to the addition of NLSls. 

test cases and concluded that uncertainties on the line fluxes we 
measured can nominally be taken as 30%. Finally, we included 
only positive measurements of line intensities, excluding a few 
cases of upper flux limits, that did not alter our basic finding as 
discussed below. We were therefore left with 130 (107) positive 
measures of the N V/C IV (N V/He 11) line ratios. Luminosities 
were derived from the monochromatic flux at A1450, taking 
Hq = 75 km s"' Mpc"^ and qa = 0.5. 

Our sample is far from being statistically complete. How- 
ever, we do not know of any observational or selection effect 
that will bias the sample against strong N V lines (that are eas- 
ily detectable and measurable). In particular, our sample of UV 
(HST) sources is strongly biased in favor of radio-loud quasars 
(e.g.. Wills et al. 1995). To test this effect, we carried out the 
analyses as discussed below while excluding line ratios of the 
radio-loud population. Our results were not significantly af- 
fected by this test (see § 3). Thus, while the following corre- 
lation slopes may depend on the specific sample, the important 
result (the fact that low-luminosity sources show very strong 
N V/C IV; see § 3) is insensitive to that. 

3. RESULTS AND DISCUSSION 

3.1. New Correlations Involving N V A 1240 

Figure 2 shows flie N V/C IV (top) and N V/He n (bottom) 
line ratios as a function of luminosity. The HF93 line ratios of 
the entire sample were also binned in ranges of 0.5 in log i/Lj, 
to minimize the effect of uneven distribution in luminosity. The 
average positions of all objects in each of those bins are shown 
as large squares with error bars in Figure 2. The error bars on 
the large squares represent, for each axis, the standard deviation 
divided by the square root of the number of objects in each bin. 
For each HF93 line Tatio-uL^, log-log diagram, we performed a 
linear regression analysis and calculated the Pearson and Spear- 

^ H/3 is the most suitable emission line for distinguishing broad-lined BLRs from narrow-lined ones. Measurements of FWHM(H/3) are not available for ~ 60% of 
our sample's objects that have z ^ 0.5. Throughout this Letter, we treat those unknown IJWHMCH/J) objects as BLAGNs. 

^ IRAF (Image Reduction and Analysis Facility) is distributed by the National Optical Astronomy Observatories, which are operated by AURA, Inc., under cooperative 
agreement with the National Science Foundation. 
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TABLE 1 

LiNLAR RLCiRLSSION PARAMhlLKS I OR N \7C I\' N'LRSIS L()(i I'L,, 



Data Set Number of Pearson Spearman Slope Constant 

Code" Objects (r) (r,) (a) (b) 

"B m OD TB 0.19 ±0.02 -2.85 ±0.23 

B+N 130 0.55 0.60 0.13±0.02 -2.11±0.23 

RQQ(B) 105 0.72 0.74 0.18 ±0.02 -2.82 ±0.23 

RQQ(B+N) 114 0.56 0.60 0.13±0.02 -2.04±0.24 

B+up.lim. 137 0.67 0.69 0.19±0.02 -2.91 ±0.24 

B+N+up.lim. 146 052 0.57 0.13±0.02 -2.17±0.24 



"Data set codes are B for BLAGNs, N for MLS Is, RQQ for radio-quiet quasars, 
and up.lim. for upper limits on the line ratio. 

TABLE 2 

Linear Regression Parameters for N v/He ii versus log uLt, 



Data Set Number of Pearson Speannan Slope Constant 

Code" Objects (r) (r,) (a) (6) 

B m rr5s fmi o. 14 ±0.02 — -1.37 ±0.27 

B+N 107 0.50 0.53 0.11 ±0.02 -0.95±0.24 

RQQ(B) 83 0.61 0.64 0.14 ±0.02 -1.25 ±0.26 

RQQ(B+N) 92 0.54 0.57 0.11±0.02 -O.87±0.23 

B+up.lim. 110 0.54 0.55 0.15± 0.02 -1.45 ±0.28 

B+N+up.lim. 119 046 0.48 Oil ±0.02 -1.00 ±026 



log nu{cm~^) < 14] that completely cover the typical range at- 
tributed to the broad emission line gas in AGNs [log nH(cin~^) w 
10]. Since there is no evidence for densities as large as 10'^- 
10^^ cm~^ in NLSls, we have no reason to suspect that they 
lie outside the range covered by the H02 calculations. Accord- 
ing to the calculations, the N V/C IV line ratio is a reliable N/C 
indicator over the range of interesting physical conditions ex- 
pected in the BLR. H02 have also considered several of the 
weak intercombination lines, such as N III] A1750/O III] A 1663 
and N IV]/C IV , previously discussed by Shields (1976), and 
concluded that the first is more reliable than the second, as it is 
less sensitive to the model assumptions. Results for all relevant 
line ratios are shown in their Figure 4 and a specific application 
to the locally optimally emitting clouds (LOCs) model is shown 
in H02 Figure 5. 

The only other line ratio available to us, except for the 
N V/C IV ratio discussed above, is N IV]/C IV. We have there- 
fore investigated the N/C obtained from the two line pairs both 
observationally and theoretically. Several BLAGNs of our sam- 
ple had published flux values (or upper flux limits) of the weak 
Une N IV]. We added to this subsample measurements of the 
N IV] line in eight of our nine extreme NLSls. The N/C abun- 
dance obtained from the N IV]/C IV line ratio was then calcu- 
lated for the subsample, assuming that upper hmits represent 
real ratios. We find that the two line ratios are weU correlated 
(r = 0.8 for 30 sources); i.e., if N V/C rv is a good N/C indica- 
tor, so is N IV]/C IV. However, the derived N/C, assuming both 
ratios are reliable metallicity indicators, is very different. The 
N V/C IV ratio gives systematically larger N/C, sometimes by a 
factor as large as 3 or 4. We stress that the results are still ten- 
tative because of the large number of upper limits rather than 
real line ratios used in the analysis. Better data are required to 
confirm this correlation. 

Regarding the suitability of N IV]/C IV as an N/C indicator, 
we note that the calculations presented in Figure 4 of H02 show 
that the line ratio does not change by more than a factor of 2 
over the range of conditions thought to be acceptable in AGN 
BLRs (ionization parameter of 0.03-0.3 for the H02 continuum 
and density below about 10'^ cm~^). The H02 conclusion that 
the line ratio is not a robust N/C indicator is based on regions in 
parameter space that are different from the one specified above. 
Moreover, the particular example shown in Figure 5 of H02, 
applicable to the LOC model, clearly shows that under such 
conditions (that produce well most of the observed fine ratios 
in AGNs; see Baldwin et al. 1995), flie N IV]/C IV fine ratio 
is indeed a very good N/C indicator A key issue is whether or 
not the physical conditions in the BLR of NLSls are similar to 
those in broader line AGNs. The idea that density and optical 
depths may be different has been proposed in the past, and a 
better assessment of the N IV]/C IV suitability in this case must 
await a more detailed theoretical investigation of such sources. 
At present we do not have a large enough sample and good 
enough calculations to test the suggestion that the results shown 
in Figure 2 are due to N V/C IV being an inadequate N/C indi- 
cator. 

3.2.2. High-Metallicity NLSls? 

The second scenario is based on the assumption that the 
N V/C IV is a reliable N/C indicator. This led HF93 to suggest a 
strong Z-L relationship in AGNs. However, our new measure- 
ments clearly show that NLS 1 s do not follow this Z—L relation- 
ship. NLSl metaUicities, as indicated by the N V/C IV ratio, are 



"Data set codes are identical to those in Table 1 . 

man Unear correlation coefficients. This analysis was carried 
out once for the BLAGN population and a second time for the 
entire data set. The results are presented in Tables 1 and 2; one 
can see that the exclusion of radio-loud objects or the addition 
of flux ratio upper limits (that were treated as real ratios) had 
Uttle effect on the correlations and slopes. 

Inspection of Figure 2 and Tables 1 and 2 shows that the 
N V/C IV ratio is strongly correlated with luminosity in the 
BLAGN case. However, it is also apparent that this relationship 
breaks down completely at low luminosities, when the NLSls 
are introduced into the sample. Several NLSls that belong in 
the lowest luminosity regime even have N V/C IV ratios as high 
as those of the most luminous high-z quasars. The extremely 
strong N V A 1240 in some NLSls has been noted earlier by 
Wills et al. (1999) who did not investigate the resulting N/C 
abundance ratio. The N V/He II ratio behaves somewhat differ- 
ently when NLS 1 s are added (Fig. 2). We find that the N V/He II 
ratio in BLAGNs is correlated with luminosity, although not as 
strong as in the N V/C IV case, and that NLSls only slightly 
deviate from the N V/He n Z—L slope. We also discover that 
both line ratios are not correlated with luminosity for vL,^ at 
1450A ^ lO"**" erg s"', neither for BLAGNs nor for the entire 
sample. We emphasize once more that the N V/C rv ratio in 
NLSls is very similar to the one observed in high-L AGNs and 
this result, which is based on an easy-to-measure line ratio, is 
enough to completely change the original HF93 correlation, re- 
gardless of the exact slope or value of the correlation coeffi- 
cient. 

3.2. High N V/C IV at Low Luminosity 

Our new results point at two distinct scenarios: either (1) the 
HF93 line ratios overpredict N/C at least under some physical 
conditions or (2) high metaUicities at low luminosities are pos- 
sible and are seen in NLS Is. We discuss briefly the implications 
of these two scenarios and defer the more detailed analysis to a 
later publication. 

3.2.1. Is the N V/C IV Line Ratio a Reliable N/C Indicator? 

Hamann et al. (2002, hereafter H02) have used state-of-the- 
art photoionization calculations to investigate several line ra- 
tios in order to select those that are robust abundance indi- 
cators. The calculations span a vast range of densities [7 < 
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similar to those of the most luminous high-z quasars in our sam- 
ple and are higher, by almost an order of magnitude, than those 
of BLAGNs with similar luminosities. In the N V/He n case, 
NLSls show only slightly higher metallicities for a given lumi- 
nosity compared with BLAGNs. This effect may be attributed 
to the more complex dependence of N V/He II on other physical 
parameters, such as the spectral energy distribution. Accepting 
this scenario, the HF93 Z-L relationship cannot be a simple 
two-parameter dependence for all AGNs. 

We also attempted to find a link between BLAGNs and 
NLS 1 s in order to see whether the Z-L relationship is a smooth 
function of FWHM(H/3). We found no correlation between 
FWHM(H/3) and metallicity, luminosity, or a combination of 
the two. Since we have FWHM(H/3) values for about half of 
our sample, we checked whether FWHM(C IV) could be more 
suitable, since this emission line is more dominant in our case. 
Again we found no correlation with the other parameters, nor 
any FWHM(C IV)-FWHM(H/3) correlation, for objects hav- 
ing both Unes measured. In fact, we find that only AGNs that 
have FWHM(H/3) < 1500 km s"' follow a significantly differ- 
ent Z-L relationship. 

If high metallicity is indeed another NLSl extreme prop- 
erty, the question whether this is related to some fundamen- 
tal NLSl physical parameter remains unanswered. The intro- 
duction of NLSls to the N V/C IV Z-L diagram completely 
changes the correlation at low L and implies that there is an ad- 
ditional dimension to this dependence that allows high metal- 
licities at low luminosities. This hidden variable in the Z-L 
relationship may be related to fundamental physical proper- 
ties, such as the accretion rate or the age of the central BH. 
The key to answering this question possibly lies in the claim 
that NLSls have low BH masses for their luminosities (equiv- 
alent to larger L/L-Edd or higher accretion rate) and therefore 
follow a different mass-luminosity relationship than BLAGNs 
(Peterson et al. 2000). Boroson & Green (1992) pointed out 
that L/LEdd might be the underlying fundamental physical prop- 
erty of their primary eigenvector. Since NLSls lie at one ex- 
treme end of this eigenvector, a natural hypothesis is that this 
property is also driving their unusually high metallicities. In 
order to test this hypothesis, one needs accurate BH mass de- 
terminations for our sample. Unfortunately, those are available 



for only 34 objects (Kaspi et al. 2000), including only three 
extreme NLSls. Moreover, measuring BH masses for high-z 
quasars is not a straightforward task since reUable determina- 
tions rely on reverberation mapping studies that would require 
at least a decade-long monitoring campaign due to cosmologi- 
cal time dilation. In addition, measurements of [O III], optical 
Fe n, and H/J lines in high-z quasars are scarce (e.g. Mcintosh 
etal. 1999). 

3.3. The Z—L Relationship in BLAGNs 

Finally, we remain within the framework of the second sce- 
nario, where the HF93 fine ratios are considered reUable abun- 
dance indicators and NLSls show systematically higher metal- 
licities than BLAGNs with comparable luminosities. We there- 
fore removed all NLSls from the sample and checked the 
Z-L relationship for BLAGNs |we caution the reader that the 
BLAGN group includes objects with unknown FWHM(H/3); 
see § 2]. Our results confirm the observational evidence for 
a correlation between metallicity and luminosity (HF93) for a 
large (statistically incomplete) sample of BLAGNs. This corre- 
lation may be written in the form Z cx L" , where Z and L are the 
BLR metallicity and BLAGN luminosity, respectively. Since 
both N V/C rv and N V/He n are approximately proportional 
to Z (HF93), the index is assigned with the mean value obtained 
for our sample, which is a ^ 0.2 (Tables 1 and 2). Combining 
this result with the empirical AGN BH mass-luminosity rela- 
tionship (Kaspi et al. 2000), Mbh oc L^-\ we find Z ex Mjin*'^. 
Alternatively, one may argue that the (Kaspi et al. 2000) re- 
sult regarding the L versus M relation is biased because of 
the small size of the sample, and a more reaUstic form is per- 
haps Mbh oc L. In this case, Z oc Mbh"^. These results dif- 
fer from those that were reached by supermassive BH growth 
considerations where a was assumed to take a value of « 0.5 
(Wang 2001). We note, again, that the above obtained a de- 
pends on the sample selection and is, therefore, rather uncer- 
tain. 

This research is supported in part by a grant from the Israel 
Science Foundation. We thank an anonymous referee for useful 
connments and suggestions. 
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